Astronomy & Astrophysics manuscript no. 


February 1, 2008 


(DOI: will be inserted by hand later) 





Pristine CNO Abundances from Magellanic Cloud B Stars 
I. The LMC Cluster NGC2004 with UVES* 

A.J. Korn 1 , S.C. Keller 2 , A. Kaufer 3 , N. Langer 4 , N. Przybilla 1 , O. Stahl 5 , and B. Wolf 5 

1 Univeritats-Sternwarte Miinchen (USM), Scheinerstr. 1, 81679 Miinchen, Germany 
e-mail: ajkorn@usm.uni-muenchen.de, nob@usm.uni-muenchen.de 

2 Lawrence Livermore National Laboratory, 7000 East Ave., Livermore, CA 94550, U.S.A. 
e-mail: skeller@igpp.ucllnl.org 

3 European Southern Observatory, Alonso de Cordova 3107, Santiago, Chile e-mail: akaufer@eso.org 

4 Astronomical Institute, Utrecht University, NL-3584 CC, Utrecht, The Netherlands 
e-mail: N.Langer@astro.uu.nl 

5 Landessternwarte Heidelberg (LSW), Konigstuhl, 69115 Heidelberg, Germany 
e-mail: ostahl@lsw.uni-heidelberg.de, bwolf@lsw.uni-heidelberg.de 



Received 2001; accepted 2001 



Abstract. We present chemical abundances for four main sequence B stars in the young cluster NGC 2004 in the 
Large Magellanic Cloud (LMC). Apart from Hn regions, unevolved OB-type stars are currently the only accessible 
source of present-day CNO abundances for the MCs not altered by stellar evolution. Using UVES on the VLT, 
we obtained spectra of sufficient resolution (R = 20 000) and signal-to-noise (S/N > 100) to derive abundances for 
a variety of elements (He, C, N, O, Mg and Si) with NLTE line formation. 

This study doubles the number of main sequence B stars in the LMC with detailed chemical abundances. More 
importantly and in contrast to previous studies, we find no CNO abundance anomalies brought on by e.g. binary 
interaction or rotational mixing. Thus, this is the first time that abundances from H n regions in the LMC can 
sensibly be cross-checked against those from B stars by excluding evolutionary effects. We confirm the H n-region 
CNO abundances to within the errors, in particular the extraordinarily low nitrogen abundance of e(N) ~ 7.0. 
Taken at face value, the nebular carbon abundance is 0.16 dex below the B-star value which could be interpreted 
in terms of interstellar dust depletion. Oxygen abundances from the two sources agree to within 0.03 dex. 
In comparison with the Galactic thin disk at MC metallicities, the Magellanic Clouds are clearly nitrogen-poor 
environments. 



Key words. Stars: abundances - Stars: atmospheres - Stars: early type - Magellanic Clouds - Clusters: individual: 
NGC 2004 



1. Introduction 

Young populous clusters like NGC 2004 in the Large 
Magellanic Cloud (LMC) are ideal objects to study both 
the current state of chemical evolution of the LMC and 
the physics of hot stars at a metallicity a factor of two 
below that of the Sun. 

The study of stellar absorption spectra yields chemical 
abundances with an accuracy comparable to that achieved 
from the analysis of nebulae such as Hn regions. However, 
both sources are subject to physical processes which alter 
the abundances of certain chemical species in a system- 
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atic fashion. In most cases, the direction of these system- 
atic effects is known, yet their extent is often uncertain. 
In particular, the scaling of these effects with metallicity 
is poorly known. In Hn regions, it is mostly dust forma- 
tion and depletion onto grains which will result in lower 
limits on the present-day (pristine) abundances. As for B 
stars, rapid rotation can lead to contamination of the at- 
mospheres with CN-cycled material from the core result- 
ing in He and N enrichment and a corresponding C (and 
- to a lesser extent - O) depletion (Fliegner et al. 1996, 
Meynet & Maeder |2000| , Heger & Langer |2000"| ) . The onset 
of "rotational mixing" is a function of the rotation rate, 
stellar mass and metallicity. Rotating models of hot stars 
indicate that it can take place very early on, i.e. whilst 
on the main sequence (MS). Furthermore, one can expect 
rotational mixing to be more efficient in metal-poor envi- 
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Table 1. Chemical abundances of bona- fide MS B 
stars in the LMC in comparison with Hii-region data. 
e(X) := log(nx/nn)+12. Typical errors are 0.2-0.3 dex. 



star 


PS 34-16 a 


LH 104-24 a 


NGC 1818/Dl b 


Hn c 


e(C) 


7.10 


7.50 


7.83 


7.90 


e(N) 


7.50 


7.70 


7.59 d 


6.90 


e(O) 


8.40 


8.50 


8.46 


8.40 


£ (Mg) 


7.00 


7.40 


7.35 




e(Si) 


7.00 


7.40 


7.10 


6.70 


e(Fe) 


7.20 




7.34 





from Rolleston etal. ( [l996| ), b from Korn et al. ( |2000| ) 
from Garnett (1999), revised in this publication 



ronments (cf. Maeder & Meynet 2001 ), as metal-deficient 
stars have smaller radii for a given mass and thus rotate 
faster for a given angular momentum. 
Among a sample of B stars only the fastest rotators 
(urot > 200 km/s) are expected to dis play c ontaminated 
CNO a bund ances (Heger & Langer 200C, Meynet & 
Maeder |2000|) . If one is interested in the pristine present- 
day abundance pattern of these elements, it is therefore 
sensible to look for and analyse the least evolved objects 
accessible. 

So far, the faintness of unevolved MC B stars (my > 15 m ) 
has prevented an extensive confrontation of the H n-region 
CNO data with that from B dwarfs. Four stars have been 



studied in detail, three by Rolleston et al. (1996) and one 
by Korn et al. (200C , hereafter Paper I). Of the three stars 
of Rolleston et al. ( 1996[ ) one (PS 34-144) is classified as 
a He-weak star which we will therefore disregard in what 
follows. The CNO abundances of the remaining three MS 
stars are confronted with the most recent measurements 
from H II regions in Table |l|: While the consistency among 
determinations of O abundances is very high, there is a 
significant positive offset (~ 0.7 dex) of B-star N abun- 
dances with respect to the H n-region value. 
It is a priori unclear which of the two data sets reflects 
the pristine LMC N abundance. There are many possible 
explanations for the peculiar behaviour of N: either the 
Hii-region value is systematically too low (note that it 
is - like all the other H n-region data - exclusively based 
on 30 Dor and N44C, cf. Garnett |l999j ), the B-star value 
systematically too high or we see the imprint of a phys- 
ical process which alters one of the two data sets. Since 
dust depletion of nitrogen seems hard to envision (Mathis 



199E), it looks as if the B stars are to blame. And indeed, 
it might be that - due to the limitations imposed on us 
by the use of 4m-class telescopes up to 1999 - choosing 
the brightest targets compatible with MS colours resulted 
in preselecting stars that are preferentially rapid rotators 
as rapidly rotating stars will evolve to higher luminosities 



(Fliegner et al. |1996| ). 

Whatever the scenario, we regard it as worthwhile to push 
the limit further towards the zero-age MS to see whether 
the N abundances already found are representative of 
LMC B stars in general. 



In contrast to the paucity of data on dwarfs, many MC su- 
pergiants - from spectral t ype O to K - have been studied 
in det ail (e. g. Haser et al 1998 : O-type , Pap er I: B-type, 
Venn |l999b|: A-type, Andrievsky et al. |200l| : F-type, Hill 
& Spite 1999: K-type). Obviously, in trying to understand 
the CNO pattern of these evolved objects the problem 
lies in separating all contributing signatures: abundance 
changes due to dredge-up episodes in the red-giant phase 
and potentially rotational mixing in earlier phases when 
the rotation rate was still high. This is particularly true 
for stars located in the so-called Blue Hertzsprung Gap 
(between the MS and the region of blue-loop excursions, 
cf. Fitzpatrick & Garmany 1990Q , whose evolutionary sta- 
tus is still highly uncertain. 

Like in the MS B stars already discussed, the most sig- 
nificant signature is th at of n itrogen: Among SMC A- 
type supergiants Venn ( 1999b ) found a large scatter in 
NLTE (non local thermodynamic equilibrium) N abun- 
dances (6.9 < e(N) < 7.7) which is hard to explain in the 
framework of the first dredge-up alone and argues in 
favour of rotational mixing in earlier evolutionary phases 
as the most likely cause. As far as abundances from 
other elements with less pronounced signatures are con- 
cerned, one has to worry about how severely they might 
be affected by systematic effects when comparing differ- 
ent spectral types: not in all cases the assumptions made 
(plane-parallel geometry, stationarity, LTE) are as justi- 
fied as in the case of dwarfs. For example, the carbon 
abundance from K supergiants seems to be systematically 
higher than in the hotter stars. On the whole, the general 
picture drawn above is, however, supported by all spectral 
types. 

All in all, rotational mixing could help to resolve the dis- 
crepancies outlined above: adding rapid rotation, a frac- 
tion of any young stellar population will display modified 
CNO abundances, widen the MS, populate regions of the 
HRD which are traversed quickly by non-rotating stars 
and lead to enhanced scatter in N after the first dredge- 
up by mimicking a distribution of initial abundances. 



2. Preselection of suitable candidates 

Since B stars are frequently fast rotators and/or show 
em ission in their Balmer lines (Be stars, cf. Jaschek et 
al. 1981 ) , a great effort was made to select the small num- 
ber of stars which are suitable for an abundance fine anal- 
ysis: non-Be and with projected rotational velocities v sini 
below 100 km/s. 

Preselection of slow rotators in NGC 2004 was 
made from spectra obtained with the Siding Springs 
Observatory (SSO) 2.3m telescope using the Double- 
beam Spectrograph. Aperture plates were fabricated 
to obtain medium resolution (0.3 A/pixel) spectra 
(4100 A < A < 5100 A) for 34 MS B stars. 



To determine v sinz we have used a x minimization 
technique which locates an optimal match between the 
observed spectrum and rotationally broadened synthetic 
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Table 2. Observation log for the programme stars. 
Nomenclature and visual magnitudes from Robertson 
(1974). i cxp denotes the expsoure time, m/y the month 



1500 



and year of the observation. The line-free wavelength 
ranges 4202-4210, 4402-4409, 4493-4500 and 4680- 
4695 A were used to estimate the S/N ratios (la) of the 
rebinned spectra. Appended to this table are the two lu- 
minosity class in stars in NGC 2004 from Paper I observed 
with CASPEC. 



object in 


mv 


C exp 


S/N 


m/y 


vsini 


NGC 2004 


[mag] 


w 






[km/s] 


B18 


14.8 


1.0 + 1.0 


140 


12/00 


130 


C8 


14.9 


1.0 


110 


11/00 


100 


C9 


15.8 


1.4 + 1.4 


100 


11/00 


60 


C16 


14.7 


1.0 


130 


11/00 


60 


D3 


15.8 


1.4 + 1.4 


120 


12/00 


70 


D15 


15.1 


1.5 


120 


12/00 


45 


B15 


14.2 


3.0 + 3.0 


70 


12/89 


25 


B30 


13.8 


3 x 2.0 


150 


11/87 


30 



NGC2004/C9 (2) (offset +500) 




4700 4800 4900 

wav elength [A] 

Fig. 1. The red part of the Blue437 exposures of 
NGC2004/C9. The blaze residuals are very constant from 
order to order allowing for a thorough rectification using 
a low-order polynomial. 



spectra interpolated from a grid of models which spans 
4000 K < T cff < 50 000 Kand 5.0 > log 5 >0.0 (Gonzales 
Delgado & L eithe rer 1996). More details can be found in 
Keller et al. ( [2002 ). Principle limitations on this technique 
are given by the S/N and resolution obtainable with the 
instrument used. 

As can be seen from Table [j], our preselection was very 
successful: four of the six targets (nearly 70 %) rotate at 
a rate low enough to analyse their metal content spectro- 
scopically. 

3. UVES observations 

Due to its high efficiency UVES (Ultraviolet- Visual 
Echelle Spectrograph) on ESO-VLT (Very Large 
Telescope) UT2 (Kueyen) is the ideal instrument for 
studying MC MS B stars (m v > 15 m ). We initially consid- 
ered using one of the standard templates for full optical 
coverage (DIC1 390 + 564nm), but rejected this option 
as it would have resulted in the loss of the strategic Si in 
triplet near 4560 A (dichroic gap). We finally chose the 
UVES standard setting Blue437nm (3760 A < A < 4980 A) 
covering all the strategic lines for a classical abundance 
fine analysis of B stars. The atmospheric dispersion 
corrector (ADC) was used and the slit oriented in such 
a way as to minimize contributions from neighbouring 
stars. The slit width was set to 1.2" (corresponding to 
a resolving power of R = 40 000) not to lose light under 
moderate seeing conditions. 

One to two spectra of each of the six preselected objects 
were exposed in Service Mode between November and 
December 2000. We make use of the pipeline- reduced 
spectra: they turn out to be of sufficient quality in terms 
of rectification, even for the use of H/3 as a gravity indica- 
tor. We note, though, that the blue part of the exposure 
(3800-4600 A) is much less well-defined in terms of the 
continuum. After rebinning to a 2-pixel resolution of 



R = 20 000 the spectra were normalized interactively and 
- in the case of two exposures - coadded. Table [t] gives a 
brief observation log including measured signal-to-noise 
ratios (S/N) and projected rotational velocities (vsini). 
Both sets of quantities fall in the range expected from 
the UVES exposure time calculator on the one hand and 
the preselection on the other. Figure [t] shows the pipeline 
reduction result for the two spectra of C9 in the region 
surrounding H/3. 

Except for B18, the heliocentric radial velocities i£f] are 
fully compatible with that deduced for the cluster by 
Freeman et al. ( |l983| ). B18 (just like B15 from Paper 
I) has a radial velocity some 40 km/s higher than the 
cluster average. These two stars do probably not belong 
to the cluster. Unfortunately, the cluster membership of 
the other programme stars is not certain either, since the 
difference in radial velocity between the cluster and the 
surrounding field is not significant within the errors. 



4. Analysis 

We utilize the three-step analysis devel oped b y and exten- 
sively discussed in Gummersbach et al. ( 1998 ) and Paper I. 
In short, log g is derived from the profile of H/3, T c s by 
balancing two adjacent ionization stages of silicon as a 
function of e(Si) and £, with £ being concordantly deter- 
mined between silicon and a set of oxygen lines. In all 
three steps, the metallicity of the underlying atmosphere 
(and with it the strength of the background opacities for 
the NLTE calculation) is varied in a self-consistent man- 
ner. This was shown to have a non-negligible effect on 
the abundances to be derived. In the simultaneous de- 
termination of T e s, log g, £, e(Si) and s(0) the average 
undcrabundance of Si and O is taken to be indicative of 
the overall metallicity [m/H]. Due to the [a/Fe] ratio in- 



herent to the MCs (see Sect. 5.2), this procedure results 
in slightly underestimated atmospheric metallicities. The 
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residual discrepancy between [m/H] and [Fe/H] is, how- 
ever, only of the order of 0.1 dex. Thus this approach is 
clearly superior to assuming [m/Fe] to be solar, as it is 
still often found in the literature. 

With regard to the NLTE line formation, new input 
physics was implemented in the model atoms for C, N and 
Mg. Here we only give a brief account of the expected reli- 
ability and refer the reader to the publications of Przybilla 
k Butler (2001) and Przybilla et al. (2001a,b). 
We refer the reader to Paper I for an error analysis and 
references to the original publication of the other model 
atoms. 

4.1. Carbon 

With equivalent widths of 160 ± 10 mA the strong C II 
feature at 4267A (multiplet 6, cf. Fig. ||) is clearly visible 
in the spectra of our MS B stars (we don't have access 
to the two strong lines (multiplet 2) in the blue wing of 
Ha). While it became clear relatively early on (Hardorp 
& Scholz 197C ) that this line yields artificially low car- 
bon abundances in LTE, a reliab le NL TE modelling was 
only achieved in the 1990s (Sigut 1996 ). Our model atom 
closely resembles that of Sigut and we follow Sigut 's ar- 
gumentation showing that the abundances derived with it 
are reliable to typically 0.2 dex. 
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Fig. 2. The loci of the programme stars with respect 
to evolutionary tracks of non-rotating stars for Z = 0.008 
(Schaerer et al. 



1993) 



Also shown are the two luminos- 
ity class in objects NGC2004/B30 and B15 from Paper 
I. No evolution of chemical elements is detected in com- 
paring the MS stars, whereas B30 from Paper I is clearly 
enriched in nitrogen and depleted in carbon. 



Departures from LTE are found to be ~ —0.1 dex, the ab- 
solute abundances have not changed in going from the old 
to the new model atom. 



4.2. Nitrogen 

The model atom has been extended to be applicable to 
early B-type stars. It includes the lowest 3 levels of N I, 
77 of Nil, the lowest 11 levels of Nm and the ground 
state of Niv. Photoionization cross-sections and oscil- 
latior strengths in Nm are adopted from the Opacity 
Project (Fernley, Hibbert, Kingston & Seaton; available 
only from the Topbas e database at http://cdsweb.u- 
5trasbg.fr/topbasc.htmj ) and electron collisons from R- 
matrix calculations by Stafford et al. ( 1994 ). 



The spectrum of nitrogen lines is far less pronounced than 
that of carbon. In fact, the strongest line at 3995 A has 
equivalent widths below 50 mA in our MS programme 
stars (cf. Fig. ||). In combination with projected rotation 
rates up to 70 km/s a clear detection requires S/N ratios of 
around 100. The situation is less critical in stars which are 
more evolved and/or whose nitrogen abundance is higher, 
e.g. NGC 2004/B30 (125 mA) or NGC 1818/D1 (74 mA, 
cf. Fig. @). 

As far as the modelling is concerned, we note that the 
subordinate lines at 4621 and 4630 A yield abundances 
fully compatible with those derived from Nil 3995. As in 
the case of carbon, these abundances are expected to be 
accurate to 0.2 dex. 



4.3. Magnesium 

The only Mg line we have access to in the B-star temper- 
ature regime is the Mgn doublet at 4481 A (cf. Fig. 0). 



5. Results 

Table | gives the stellar parameters and abundances de- 
rived for the four programme stars in comparison with 
H II- region data from Garnett ( 1999| ) . Consistency among 
the abundances is generally very high with the rms scatter 
confined to well below 0.1 dex. Also included is the main 
sequence star Dl in the LMC cluster NGC 1818 already 
analysed in Paper I. For the latter star, the abundances 
of He, C and N are redetermined on the basis of all avail- 
able lines and utilizing the new model atoms. Where avail- 
able, the abundances are compared with measurements 
f rom L MC H n regions as carried out / compiled by Garnett 
( |1999| ). Finally, we also in clude the Solar System compo- 
sition (Grevesse & Sauval 1998 ) as a reference. Note that 
a recent redetermination of the solar oxygen abundance 
from [O i] 6300 based on 3D h ydrodynamical computa- 
tions (Alle nde P rieto et al. |200lD yields e(O) = 8.69 ± 0.05. 
Holwcgcr ( 2001 ) comfirms this trend towards a lower value 
by deriving e(O) = 8.74 ±0.08 from eight optical and IR 
lines taking into account effects due to NLTE and granula- 
tion. Thus the difference between oxygen abundances from 
stars in Orion and the Sun reduces to 0.1 dex at most (cf. 
Gummersbach et al. 1998| ), but important consequences 
can also be expected for the evaluation of the [a/Fe] ra- 
tios in the MCs (see below). 

Figure |^ shows the positions of the programme stars in a 
log T c ff - log g parameter space, sometimes referred to as 
Kiel diagram. Note that the four MS stars nicely span 
the full width of the main sequence in the B-star temper- 
ature regime. This fact allows us to potentially address 



A.J. Korn et al.: Pristine CNO Abundances from Magellanic Cloud B Stars 



5 



Table 3. Derived stellar parameters and abundances for the four programme stars plus NGC1818/D1 from Paper 
I. Solar System composition from Grevesse & Sauval (1998), LMC Hn-region abundances from Garnett ( 1999] ). The 
customary notation e(X) := log(nx/^H) + 12 is employed. Representative error estimates are assigned to the individual 
abundances determinations. 



object 




logs' 




VSVQ.I 


e(He) 


e(C) 


e(N) 


e(O) 


E(Mg) 


e(Si) 


e(Fe) 




[K] 




[km/s] 


[km/s] 














LTE 













10.99 


8.52 


7.92 


8.83 




7.55 


7.50 


NGC2004/D15 


22 500 


3.80 





45 


10.95 


8.04 


6.95 


8.29 


7.35 


7.04 


7.30 


± 


1000 


0.20 


2 


10 


0.2 


0.2 


0.2 


as 


as 


as 


a 5 


NGC2004 /C16 


24600 


3.95 


4 


60 


10.94 


8.10 


7.05 


8.43 


7.29 


7.10 


7.30 


± 


1000 


0.20 


S 


10 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


a 5 


NGC2004/D3 


23 900 


4.15 


1 


70 


11.00 


8.15 


7.05 


8.40 


7.43 


7.20 


7.35 


± 


1000 


0.20 


2 


10 


0.2 


0.2 


0.2 


as 


as 


as 


0.3 


NGC2004/C9 


24 000 


4.35 





60 


10.93 


7.95 


7.00 


8.34 


7.40 


7.05 


7.35 


± 


1000 


0.20 


2 


10 


0.2 


0.2 


0.2 




0.2 


as 


a 5 


^(-^-)2004 










10.96 


8.06 


7.01 


8.37 


7.37 


7.10 


7.33 












0.03 


0.09 


0.05 


0.06 


0.06 


0.07 


0.03 


H II LMC 










10.95 


7.90 


6.90 


8.40 




6.70 




± 










0.2 


0.2 


0.2 


as 




as 




NGC1818/D1 


24 700 


4.00 





30 


11.04 


7.80 


7.40 


8.46 


7.35 


7.10 


7.37 


± 


1000 


0.20 


2 


10 


0.2 


0.2 


0.2 


as 


as 


as 


0.3 



chemical evolution of certain elements in these stars dur- 
ing their MS lifetime: no such evolution is, however, found 
(cf. Table | in which the stars are ordered according to 
their gravity). 

Figure |^ gives examples of profile fits for the most impor- 
tant species of this study: H, CNO, Si and Mg. NLTE 
effects in iron are currently under investigation. Thus, 
for the time being, e(Fe) based on the few Fein lines 
(Wa < 40 mA) visible between 4000 A and 5000 A is given 
assuming LTE. 



5.1. Pristine CNO abundances in the LMC 

In comparing the stellar CNO abundances with the neb- 
ular ones, one immediately notices several features: the 
oxygen abundances are in excellent agreement, while both 
nitrogen and carbon are somewhat more abundant in the 
stars. As carbon is one of the principle dust-forming ele- 
ments, interstellar dust depletion is the first reason that 
comes to mind. Though not of high statistical significance, 
this offset of 0.16 dex could be interpreted as a manifes- 
tation of the fraction of carbon locked up in interstellar 
dust. Unfortunately, the error limits on both sides prevent 
us from drawing firm conclusions on this matter. Another 
observation, however, may lend support to this scenario: 
a very similar offset is found between the gas-phase car- 
bon abundance of Orion (—8.4, cf. Garnett et al. 1995) 
and that of the Sun (between 8.52 (Grevesse & Sauval 
199S) and 8.59 according to Holweger 2001). More work 



Conversely, nitrogen is believed not to de plete into the 
dust phase in significant amounts (Mathis |l 996 ) . An en- 
hanced nitrogen abundance in B stars could then be taken 
as an indication of mixing (internal or external), raising 
the original (nebular) abundance from which the stars 
formed some million years ago. The high internal consis- 
tency among the e(N) determinations in the programme 
stars speaks against this scenario as a likely explanation. 
Besides, the residual offset with respect to the nebular 
abundance is hardly of high significance considering the 
accuracy of 0.2 dex for each of the data sets. We there- 
fore take this result as strong evidence in favour of the 
extraordinarily low present-day nitrogen abundance 
of e(N) ~ 7.0. On the square-bracket scalef] this corre- 
sponds to [N/H] < -0.92 at a metallicity of [Fe/H] ~ -0.2 
or [N/Fe] ~ -0.7. [N/Fe] ~ is only reached after the first 
dredge-up on the red-giant branch, as already emphasized 
by Venn ( |1999a| ). 

Finally, we comment on oxygen-related matters: there 
can be little doubt that the present-day LMC oxy- 
gen abundance is close to 8.40. With a mean nitro- 
gen abundance of 7.01, our mean oxygen value of 8.37 
results in log(N/0) = — 1.36. Therefore, even though the 
metallicity of the LMC is not so different from solar, its 
N/O ratio is among the lowest found in the local universe, 
as derived from Galactic and extragalactic nebular and 
stellar abundances (cf. Henry et al. 2000 ). Indeed, within 
the chemical evolution model favoured by Henry et al., 
which attributes nitrogen production entirely to low and 



is clearly needed here. 



[X/Y]*:=log(2t) 



l°g(^)e 
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intermedate-mass stars, this implies that the present-day 
nitrogen abundance of the LMC is still dominated by pri- 
marily produced nitrogen, which is consistent with the 
metallicity-dependent nitrogen yields of van den Hoek & 
Groenewegen ( |1997| ). 

5.2. The a-elements 

Due to a state-of-the-art redetermination of the solar oxy- 
gen abundance (see above) it has become unclear what 
[O/Fe] ratio s(0) = 8.40 corresponds to. In fact, this ques- 
tion concerns not only oxygen, but all other a-elements 
(Ne, Mg, Si, S, Ca and possibly Ti) in the sense that oxy- 
gen is often taken to be indicative of the behaviour of these 
elements. Thus the whole question of whether or not the 
[a/Fe] ratios in the MCs are vastly different from those 
found in the Galactic thin disk at the appropriate metal- 
licities must be re-evaluated in the light of these recent de- 
velopments: From the three a-elements we have access to, 
we derive [a/Fe] ~ —0.15 using the new solar oxygen abun- 
dance. Even if we account for a lower NLTE Fe abundance, 
[a/Fe] is still significantly lower than in t he G a lactic thin 
disk: +0.10 at [Fe/H] ~ -0.3 (Fuhrmann |1998| , |2000[ ). 
For silicon we find a relatively clear offset of stellar silicon 
abundances relative to that from Hll regions of 0.4 dex 
which could be significant. The LMC giants of Paper I 
show a similar offset (~ 0.3 dex) which can most naturally 
be explained as the fraction of silicon locked up in inter- 
stellar grains (silicates). 

5.3. Iron 

On average, the iron abundances presented in Table |^ are 
merely 0.17 dex below the solar value which is quite a high 
value for LMC objects. We reemphasize, however, that 
LTE results are given in Table |^. First results of the new 
NLTE computations indicate that the LTE results consti- 
tute upper limits to the NLTE abundances. Thus, the high 
iron abundances is most likely caused by the inappropri- 
ateness of LTE for stars as hot as 25 000 K. For the time 
being, we account for this fact by assigning larger error 
bars (±0.3 dex). 

5.4. Isochrone ages 

In comparison with evolutionary tracks, the stellar para- 
meters may serve to constrain the cluster age which has 
only been determined from colour-magnitude diagrams 
(CMD) so far. Biases are certainly present in both cases: 
MS fitting in CMDs will be prone to Be stars and rapid 
rotators "above" the MS mimicking a hotter turnoff and 
artificially lowering the derived age (cf. Grebel et al. 
1996). Additionally, (differential) reddening is always an 
issue. In the case of spectroscopic ages derived from in- 
dividual stars, unknown cluster membership might intro- 
duce a bias, not even the direction of which is known. 
Advantageously, reddening does not enter. 



Having shown that our programme stars are bona-fide un- 
mixed objects, we are hopeful that the derived ages are at 
least not affected by rotation. As can be seen from Fig. |[ 
the ages scatter substantially, e.g. C9 is fully compatible 
with age zero. It is unclear whether or not this hints at 
multiple star formation episodes for NGC 2004, as pro- 



posed by Caloi & Cassatella (1995) based on IUE data. 
Our best estimate comes from the most evolved MS tar- 
get D15 close to the turnoff where the age resolution is 
highest: 19Myr < t N GC 200 4 < 24Myr. A CMD analysis 
by Bencivenni et al. (1991) points towards a much lower 



age of 8Myr confirming an earlier result by Hodge (1983). 
Note, though, that these estimates are based on B — V 
colours which are a rather insensitive indicator of turnoff 
temperatures and prone to conta mination with Be stars. 
More recently, Keller et al. ( |2000|) utilized the HST filters 
F160BW (far UV), F555W (V) and F656N (Ha), sepa- 
rated the Be population and derived an age for the cluster 
of 16±4Myr, in much better agreement with our result. 

5.5. A fresh look at NGC1818/D1 

Having redetermined the He, C and N abundance of this 
star, it is worthwhile taking another look at it: even 
though both e(He) and e(N) have decreased due to better 
statistics and more sophisticated modelling, the signature 
indicative of mixing with CN-cycled material remains in- 
tact. More importantly, we now have pristine B star abun- 
dances to compare it to (see Fig. ||). 

The agreement in the O, Mg, Si and Fe abundance with 
NGC 2004/C16 (the object closest to Dl in terms of stellar 
parameters) is remarkable, yet both carbon and nitrogen 
are offset by > 0.3 dex in opposite directions. We refrain 
from rephrasing the discussion already presented in Paper 
I, but simply stress that NGC1818/D1 remains a prime 
candidate for a MS B stars that has undergone rotational 
mixing. We urge observations of boron to be perfo rmed 
for this star which according to Fliegner et al. ( 1996 ) have 
the potential of distinguishing between external contam- 
ination (e.g. binary interaction) and an internal process 
like rotational mixing. 

6. Conclusions 

We have presented detailed chemical abundances for four 
MS B stars in NGC 2004. This study doubles the num- 
ber of MS B stars in the LMC for which accurate abun- 
dances have been derived by means of high-resolution 
spectroscopy. None of the stars shows any indication of 
effects due to binary interaction or rotational mixing al- 
lowing us to derive what we consider to be the first data 
set of unaltered present-day CNO abundances from LMC 
stars. 

— We confirm the extraordinarily low LMC nitrogen 
abundance previously found from H 11-region studies to 
within 0.1 dex. This value is more than 0.5 dex below 
average values found in the Galactic thin disk at LMC 
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metallicities (Liang et al. 2001, cf. their Fig. 10). It im- 
plies an enrichment history for the two environments 
different from one another and - within the chemical 
evolution model of Henry et al. ( 200C ) - a dominance 
of primary nitrogen production in the LMC until 
today. 
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depletion would act, the offset found is hardly signif- 
icant considering the error limits on both sides. More 
work is needed to reduce both random and systematic 
errors to draw definitive conclusions on dust depiction 
fractions. 

— The B-star oxygen abundance (also from B giants, cf. 
Paper I) is in excellent agreement with the nebular 
value. In the absence of systematic offsets between the 
two data sets, little nebular oxygen seems to be tied 
up in grains at LMC metallicities. 

— The stellar silicon abundance is offset from the nebular 
one by 0.4 dex. This might imply that a major fraction 
of the interstellar silicon in the LMC is bound in grains. 

This work illustrates the potential of abundance analyses 
of hot stars based on high-quality spectra and sophisti- 
cated input physics. With modern multi-object spectro- 
graphs (e.g. GIRAFFE on VLT UT2 going on-line in 2002) 
we will be able to analyse some 100 slowly rotating B stars 
in each Galactic and MC cluster like NGC 2004. This will 
not only result in reduced random errors (which we will be 
able to address truly statistically for the first time), but 
also prove or refute the importance of rotational mixing 
for this class of stellar objects. In particular, a sequence 
consisting of Galactic, LMC and SMC B stars could settle 
the role metallicity, rotation and mass loss play in hot- 
star evolution. As for the Hn-region data, more objects 
will have to be observed to beat down the errors on that 
side as well. 

Systematic errors — arguably more important than ran- 
dom ones at the current level of accuracy - will have to 
be addressed by an integral approach: by studying the 
stellar and gas component of H n regions in tandem. Both 
fields of research would undoubtedly profit from such an 
endeavour. 
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Fig. 3. Profile fits for H (gravity determination), Siii/ni/iv (temperature and microturbulence indicator), Cn, Nil, 
On (microturbulence indicator) and Mgn. The observed noise level can fully account for the difference between theory 
and observation. Note that the He abundance is not based on He I 4121 A, since the observed profile seems to contain 
unidentified blend contributions other that from helium and oxygen 
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Fig. 4. Carbon and nitrogen profile fits for NGC 1818/D1: the best-fit value is compared to the mean value as derived 
from the four MS stars in NGC 2004. The thin line represents the LTE utility SYNSPEC (T c g = 25 000 K, log g = 4.0, £ = 
lOkm/s, solar composition, cf. http: / /fcros.lsw.uni-hcidclbcrg.dc/cgi-bin/wcbsynspcc.cgi[ ) which is used to locate the 
continuum. Note that some of the line strengths are greatly overestimated under LTE conditions (e.g. Cu 4616/4619 
and 4625). 



